Electron-phonon bound state in graphene 
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The fine structure of the Dirac energy spectrum in graphene induced by electron-optical phonon 
coupling is investigated in the portion of the spectrum near the phonon emission threshold. The 
derived new dispersion equation in the immediate neighborhood below the phonon threshold cor- 
responds to an electron-phonon bound state. We find that the singular vertex corrections beyond 
perturbation theory increase strongly the electron-phonon binding energy. The predicted enhance- 
ment of the effective electron-phonon coupling can be measured using angle-resolved spectroscopy. 



Discovery of graphene [TJ [2] with its unique coni- 
cal gapless bandstructure provides a new rich area for 
investigations of many-body physics of chiral massless 
fermions [3]. Active theoretical efforts are directed to 
the study of the interaction effects of Dirac carriers with 
elementary excitations of charge density waves [U [5] 
and lattice vibrations [6--8J, which result in the veloc- 
ity renormalization of the bare Dirac spectrum [9Tfl3] 
and in the formation of new quasiparticles such as plas- 
marons P3J [TS] , polarons [THl E] , and plasmon-phonon 
complexes [18 . The spectral and damping properties 
of these quasiparticles have been studied by means of 
such powerful experimental tools as Raman [30] and 
angle-resolved photoemission [21] spectroscopies. These 
high precision measurements indicate significant modifi- 
cations of the peculiar graphene bandstructure, induced 
by electron-electron and electron-phonon interactions. 

The aim of the present paper is to investigate how 
the Dirac spectrum in graphene is modified near the 
threshold of optical phonon emission. Previous stud- 
ies of this problem in bulk semiconductors |22| have 
shown that even for weak electron-phonon coupling, a, 
the branches of the true spectrum can be classified 
into three main groups, according to the effective av- 
erage number of phonon states, N = 1 — Z , bound to 
the electron and forming complex quasiparticles. Here 
Z = (1 — 9S(e,p)/9e) _1 | e=e ( p ) with E(e,p) the electron 
mass operator and s(p) the energy of the electron-phonon 
complex with momentum p. The first group includes 
states with N ~ a 2 . The spectrum of these branches 
differs from the bare one by a simple renormalization of 
the energy of the order of a 2 . The second and third 
groups are the hybrid and bound states of electrons and 
phonons. The universal threshold nonanalyticities [23] 
are responsible for the existence of these electron-phonon 
complexes, respectively, with N ~ 1/2 and N ~ 1. Note 
that for the formation of the hybrid states, a resonance 
situation in the bare spectrum is required. 

The character of the threshold singularities depends on 
the competition between the kinetic energy of the elec- 
trons and phonons and their interaction energy, there- 
fore, is largely determined by the dimensionality of the 
system. In graphene because of its two-dimensional na- 
ture the size-extend of electrons and phonons is strongly 



reduced and one can expect that the threshold nonanalyt- 
icities will result in strong modifications of the bare Dirac 
spectrum. Such a strong enhancement of the electron- 
phonon effective coupling and the binding energies of 
complex quasiparticles has been previously found in two 
dimensional semiconductor structures [2H US] . 

Recent perturbative calculations [5HTT1 ITo] of the elec- 
tron mass operator, S(e), in graphene showed that in 
rt-doped samples the real part of the lowest order contri- 
bution, So( £ )i diverges logarithmically at energies near 
the optical phonon emission threshold, e c — ±ujq, while 
the imaginary part of E ( £ ) ? related to the single-particle 
density of states, makes a discontinuous jump at the same 
threshold points. The + (— ) sign refers to the phonon 
emission process by Fermi electrons (holes), loq = 196 
meV is the longitudinal optical phonon energy. We will 
use h — 1 units. Due to this non-analyticity the effective 
electron-phonon coupling becomes large and a pertur- 
bative calculation of £o(e) in the neighborhood of the 
phonon emission threshold will not be a good approxi- 
mation for In terms of the diagrammatic expan- 
sion, the diagrams with dangerous intersections along 
one phonon and one electron lines are responsible for 
the threshold singularities [25]. The simplest diagrams 
for with such dangerous intersections are shown in 
Figs.ga-c). As seen, the number of dangerous intersec- 
tions increases with the order of the diagram and pertur- 
bation theory fails to converge when e — > e c . In order 
to find the true spectrum of the system in the neigh- 
borhood of e c , it is required to sum an infinite number 
of divergent diagrams with such dangerous intersections. 
For this we exploit an approach beyond perturbation the- 
ory [22] , which leads to coupled integral equations for the 
exact electron Green function and for the exact electron- 
phonon vertex part, drawn in Fig. [2] By solving these 
equations we find a new dispersion equation for the com- 
plex quasiparticle energy that in the immediate neighbor- 
hood below the threshold describes an electron-phonon 
bound state in graphene with new analytical dependen- 
cies on the electron-phonon coupling. The true spectrum 
does not asymptotically tend to the phonon energy but 
always remains below it at a small but finite distance. At 
finite wave vectors we find that the singular vertex cor- 
rections to So( e ) increase strongly the electron-phonon 
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FIG. 1. The simplest diagrams for the electron self energy: 
(a) the lowest order mass operator Eo(e), (b) and (c) are 
vertex corrections to Eo(e) in next two orders. The vertical 
dash-dotted lines show the dangerous intersections along one 
phonon and one electron lines corresponding to the threshold 
singularities of the mass operator, (d) The simplest diagram 
for the electron-phonon four vertex part from Fig. [2] 



binding energy in comparison with that calculated within 
the perturbative approach. 

We search for new branches of electron-phonon com- 
plex states from the poles of the exact single-electron 
Green function, G^(e), in the total energy parameter e. 
Usually in doped graphene samples the Fermi energy, £p, 
is much larger than the lattice temperature T and even 
at room temperatures most of the electrons are below £p. 
In the chiral basis the zero temperature Green function 
of noninteracting electrons 

G Je,k) = 1 (1) 

; e + EF-e^k+iO-sgne v ; 

corresponds to the thin solid lines in Figs. [T] and [2] The 
electron bare energy in the vicinity of the Dirac points in 
graphene has a linear dispersion, = /j,vp\k\, where 
k and vf are the momentum and the Fermi velocity 
of the massless fermions, described by the Hamiltonian 
Hq = —vpa-k. The Pauli matrices a act in the pseu- 
dospace of graphene sublattices and \i = ±1 labels the 
electron chirality. At low T absorption of phonons by 
electrons is negligible and we replace all exact phonon 
Green functions (dashed lines) by free phonon propaga- 
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where cj S q and q are the energy and momentum of the s 
phonon mode in graphene. 

The small dots in Figs. [T] and [2] correspond to the bare 
electron-phonon vertex functions 
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where the electron wave functions for the K point 



are VVk( r ) = (/z e l ^ k ) exp(ikr)/V-4 with A the nor- 
malizing area and 4>k the polar angle of the vector 
k. The perturbation of the graphene lattice poten- 
tial created by a single sq phonon mode is given by 
Vsq(r) = avFV s (q) exp(iqr)/v^4 with the interaction 
matrices represented as [57] 
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FIG. 2. (top) Dyson equation for the exact electron Green 
function G M (e) (thick solid lines), (bottom) The ladder type 
equation for the exact electron-phonon vertex part T^ M (e) 
(large bold dots). The solid square represents the irreducible 
four vertex part D(e) with two electron and two phonon ex- 
ternal lines. The simplest diagram corresponding to D(e) is 
drawn in Fig. |TJd) . 



for the longitudinal (s = 1) and transverse (s = 0) 
phonons. The dimensionless coupling constant is defined 
as a — (3/b 2 y / 2aui scl where a is the surface mass density 
unit cells, j3 ~ 2 a dimensionless tight-binding parame- 
ter, b the bond length between adjacent carbon atoms. 
This yields weak coupling with a 2 ~ 0.02 [TT] , 

Further we focus only on the part of the spectrum near 
the threshold of longitudinal optical phonon emission by 
electrons, e c = +oj . The singular behavior of the spec- 
trum near the threshold of phonon emission by Fermi 
holes, e c = —u>q, as well as for the transverse optical 
phonons can be treated independently in a similar way. 

In the energy region of our interest, e e c , we can 
make several simplifications using the threshold approx- 
imation. In equations corresponding to the diagrams of 
Fig. [2] we take the electron Green functions as retarded 
and after the integration over the phonon energetic pa- 
rameter oj replace it by u>q in all internal electron lines. In 
proximity of the threshold in the conductance band, the 
leading contributions to the summation over the chirality 
of the internal electron lines make the singular terms with 
p, = +1 and p! = +1. The energetic parameter of the ex- 
act Green functions in the dangerous intersections in the 
second terms in the rhs of the equations in Fig.[2j e — cj , 
lies far from the threshold e c where perturbation theory 
is applicable. Therefore, the internal exact Green func- 
tions can be replaced by the bare function Go+(e — Wo, k). 
Furthermore, in all integrations over the absolute values 
of the electron momenta corresponding to the dangerous 
intersections in Fig. [2j only the small regions near the 
Fermi wave vector, |k| ps kp and |k'| pa fcp, give the main 
contribution to the integrals in these equations. There- 
fore, in this approximation one can take the quantities 7, 
□ , and r out of the integrations over k and k' . This al- 
lows us write the Dyson equation corresponding to Fig. [2] 
as 

G-\e,p) = G^(e,p) - A+(e) f" dfa (4) 
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FIG. 3. (right) The energy spectrum of the electron-optical phonon quasiparticle for the electron density n — 5.6 x 10 cm" . 
The horizontal solid line is the phonon emission threshold, (mid) The corresponding density of states g{e) of the quasiparticle 
in units of g F = 2kp/nvF- (left) The electron-phonon binding energy versus the bare coupling g for p — pf = 1.5go- Inset 
shows the binding energy of the hybrid states at the resonance p — pf — qo- In all figures the dashed lines represent the bare 
Dirac fermions, the dotted and dot-dashed curves are calculated, respectively, within the Rayleigh-Schrodinger and Wigncr- 
Brillouin perturbative approaches. The solid curve is obtained within the present theory, taking into account the singular 
vertex corrections beyond perturbation theory. 



and represent the equation for the nonanalytical vertex 
part in the following way 

T +M (e |k F , p; p k F ) = 7+AI (k F , p; p - k F ) (5) 

+ A+(e) / |k P ,p- k F ;k F ,k F - p) 

Jo 

x r +M (e|k F ,p;p-k F ) . 
Here we introduce the following singular function 
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where the unit step 9(x) function is what remains from 
the Fermi functions at T = 0. The divergence of the 
integral at large values of k is related to the linearity of 
the graphene energy band. The cutoff of the integral at 
large momenta of the order of the inverse lattice constant 
contributes to the regular part of the integral. We are 
interested in its singular part due to the phonon emission 
threshold. It comes from the low limit of the integral, i.e. 
from the momentum range close to the Fermi wave vector 
kp, which gives 
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Because the energetic parameter of the electron inter- 
nal Green function Gq^" (e — 2wo) in the four vertex part 
□ ++ (e) is far from the threshold e c (the four vertex part 
has no dangerous intersection along one electron and one 
phonon lines), \3 ++ (e) can be expanded with respect to 
a and be replaced by the simplest diagram shown in 



Fig. [ljd). We restrict ourselves to highly doped sam- 
ples where the Fermi energy is larger than the phonon 
energy. In this regime the important contribution to the 
four vertex part D++(e) comes from scattered virtual 
phonons with approximately equal antiparallel momenta 
q = p k and q' = p k' and we replace £ p - q '- q by 
s F in the electron Green function and retain only the 
term with chirality ji" = +1 in the sum corresponding to 
the internal electron line. For such dominant scattering 
events the vertex parts 7 in D ++ (e) depend only on q or 
q' and the four vertex part becomes decoupled as 

□ ++ (w |q;q / ) « G 0+ (-o;o,fc F )7 ++ (q)7 ++ (-q') . (8) 

Taking /1 = +1 in all electron external lines and intro- 
ducing a new amplitude 
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#q7++(-q) r ++(^|p- q,p;q) (9) 







as well as the form factor 
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and with the help of Eqs. ^-§K>§ we get 

T(P) 



r(e,p) = 



l-A + (e)G a+ (-u 0l kp)r( P ) 



(10) 



(11) 



Combining Eqs. ([9|, and ( 11 ) and making use of the 
explicit expressions for the respective functions in these 
equations, we derive the following dispersion relation 
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Here we introduce a renormalized electron-phonon cou- 



pling g = a ep /(jJq. Eq. (12 1 includes the singular ver- 



tex corrections beyond perturbation theory and results 
in new properties of the electron-phonon quasiparticle in 
graphene. The perturbative Wigner-Brillouin (WB) ap- 
proach is recovered if the dominator in the rhs of Eq. ( 12 ) 



is set equal to 1. Within the Rayleigh-Schrodinger (RS) 
theory the exact energy in the rhs of Eq. ( 12 ) addition- 
ally should be replaced by e — vp(p — pp). The cor- 
responding single-particle density of states is given by 
g(e) = Ap/ (2irde / dp)\ e=e ^ where e(p) is the solution of 
Eq. @2). 

In Fig. [3] we plot the spectrum characteristics of the 
electron-phonon quasiparticle, calculated from Eq. (12 1. 
We consider only the energy range < £ < ujq . Above the 
threshold, e > uj , there exists a continuum of decaying 
states and no true elementary excitations. For e < 0, the 
threshold e c = — uo, which refers to the Fermi holes, has 
a similar effect on the spectrum. As seen in Fig. |3jleft) 
at small values of p — pp < qo (qo = ojq/vp) the effect 
of vertex corrections is weak and the spectrum obtained 
from perturbative RS and WB theories provides an ad- 
equate description. In this limit the average number of 
phonon states N ~ g and the density of states is linear 
with e. At e = the shift of 17(e) from gp, the density 
of states of the bare Dirac fermions at the Fermi level, 
is proportional to the coupling g (Fig. |3Fmid)). With 
an increase of the momentum p the RS approach fails 
completely while the spectrum, obtained within the two 
other approaches for p — pp ~ qo, describes the electron- 
phonon hybrid states with N ~ 1/2. It is seen in the 
inset in Fig. [3jleft) that the vertex corrections increase 
the binding energy, = uio — £, but the effect is still 
relatively modest. 

For even larger momenta p — pp > qo the spectrum 
obtained within WB perturbation theory (this approach 
actually has been used in Refs. I9HTT| and [16} converge 
asymptotically to the phonon energy when p — > 00 hence 
does not support an electron-phonon bound state. In 
contrast, the vertex corrections become especially im- 
portant in this region and open a small gap under the 



threshold. The true spectrum obtained from Eq. (12 1 al- 



ways remains below the phonon energy and corresponds 
to the electron-phonon bound state with N ~ 1. The ve- 
locity of the bound state tends to zero while the density 
of states increases strongly with p. Although the bind- 
ing energy of the bound state remains finite for p — > 00, 
it exhibits a stark exponential dependence on the cou- 
pling constant, = ep exp(— 1/g), and with its sub- 
Kclvin value for experimentally accessible values of g 
is hardly measurable. At finite momenta p — pp > qo, 
we find, however, that the singular vertex corrections in- 
crease strongly the binding energy e^ in comparison with 
that obtained within the perturbative WB approach (see 
Fig. [3^1eft) ) . For p — p F = l.bqo and for the doping level 
n = 5.6 x 10 13 cm~ 2 (corresponding to g sa 0.09) we find 



for the binding energy ss 14.5 meV. It is about a factor 
of 5 larger than the corresponding WB perturbative value 
and this difference increases strongly with g (Fig.[3]deft)). 
This enhancement results in a significant deviation from 
the linear Dirac spectrum that should manifest itself in 
angle-resolved measurements with the resolution smaller 
than 10 meV as stark delta-function peaks at frequencies 
ujq — and wave vectors larger than qo- 

In conclusion, we have calculated the fine structure of 
the Dirac spectrum in graphene in the proximity of the 
phonon emission threshold. The renormalized spectrum 
in the immediate neighborhood below the threshold cor- 
responds to the electron-phonon bound state. Our cal- 
culations result in a strong enhancement of the electron- 
phonon binding energy due to the singular vertex correc- 
tions, which can be probed in experiment. 
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